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Abstract 



A brief summary is given of some recent perturbative QCD results on the eval- 
uation of cross sections and on the the structure of final states in processes with 
incoming hadrons at small x. A new Monte Carlo simulation which includes 
these theoretical features is described together with its applications to the study 
of heavy quark production in ep collisions at very high energy. 



Heavy flavour production is a process with a hard scale given by the quark mass M. At very 
high energies this scale is much smaller than the collision cm. energy ^/s. This implies that 
together with the collinear logarithms (powers of lnM 2 /A 2 ) we must resum also powers of 
In a; with x ~ M 2 /s <C 1. Much progress has been made over the past few years in the 
theoretical understanding of small-x processes. In particular the key acievements in the 
region x — > are: (i) a better understanding of the "Lipatov" anomalous dimension ffl, E3]; 

(ii) the resummation || to all loop order the leading contributions in the coefficient function; 

(iii) the extension to this region of the coherent branching process; (iv) the possibility to 
resum all these new contributions by Monte Carlo methods @, |5], ||. In this talks I will 
briefly summarize these theoretical results and describe a recent application of the Monte 
Carlo simulation to heavy flavour leptoproduction at Hera and higher energies. 



1) Structure function. 

The structure function is given in term of the space-like anomalous dimension 7^ (as) (the 
limit x — > corresponds to N — > 1, where N is the energy moment index). The leading 
contributions in 7^(05) are given by an expansion in powers of cts/(N — 1) known since 
long time and recently studied in the framework of hard processes . The first terms 
of the "Lipatov" anomalous dimension are 

where as = CaolsI^ and is the Riemann zeta function. There ere are no leading terms 
of order a s , a s , and a%. Although each term is singular only at N = 1, this expansion 
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develops a square root singularity at N = 1 + (4 In 2)0:5. The presence of this singularity 
at N > 1 implies that the behaviour of the structure function for x — > is more singular 
than that given by any finite number of loops. For fixed as and small x the behaviour of 
the one loop structure function is 

xF {1) (x,Q) ~exp y/ala(l/x) , a = Aa s ln(Q 2 /Q 2 ) . (2) 

By summing the all loop result in Eq. (P one finds instead the following behaviour 

xFW(x,Q) ~aT p , p=(41n2)a s , (3) 
which is much more singular for x — > 0. 

In spite of this quite different behaviour, it turns out that at this highly inclusive level the 
all- loop and the conventional one-loop formula give similar results || . This is partially due 
to the fact that the first correction to the one-loop expression of the anomalous dimension is 
to order a s . Thus the steeper behaviour of the structure function is seen only for very low 
x. Moreover it has been pointed |S| that the steeper behaviour for as running is even more 
asymptotic than for fixed as- This is due to the presence of the cutoff in the exchanged 
transverse momenta k ti > Q . Although this condition is asymptotically negligible, it has 
some effect in reducing the evolution of the branching in the first steps. As a result the 
distribution at small x is somewhat reduced. For a detailed discussion see || and the 
contribution by Ryskin to this conference. 

2) Cross section. 

To compute the heavy flavour production cross section one convolutes the structure function 
with the coefficient function. Recently a way has been found to sum all leading terms 
ag/(N — l) n in the coefficient function as well. The key development || in this study has 
been a generalized factorization theorem in which one takes into account the off-shellness 
of the hard scattered parton. One introduces the hard elementary off-shell cross section 
<r(k t , Q) for a photon and a gluon to produce a heavy quark- ant iquark pair, where q 2 = —Q 2 
and k 2 = —k 2 are the photon and gluon squared masses respectively, and the generalized 
proton structure function J-"(x,^,k t ) giving the probability (per unit of Inx) of finding a 
gluon at longitudinal momentum fraction x and transverse momentum k t in a hard process 
at the scale //. Integrating this distribution over k t < \x one obtains the gluon structure 
function F(x,fi 2 ). By studying a(k t ,Q) at Hera energy and M = 5GeV one has that 
when W 2 3> M 2 ^> Q 2 ^> A 2 , W being the hadronic cm. energy, the natural cutoff is 
around 4M 2 , while for deep inelastic production with W 2 S> Q 2 S> M 2 the cutoff becomes 
Q 2 . Therefore the hard scale is typically assumed to be fi 2 = 4M 2 + Q 2 . However at small 
Q 2 or at W ~ M, the dynamical suppression in k t is at a significantly smaller scale. 

In the conventional calculation, the heavy flavour cross section is obtained by convolut- 
ing the on-shell elementary cross section <r(0, 0) and the gluon structure function F(x, fi 2 ). 
This procedure has two main effects: (i) for k\ < /j, 2 the elementary cross section is over- 
estimated by its on-shell value <r(0, 0); (ii) the 'tail' of the cross section at k 2 > /1 2 is 
ignored. 

Asymptotically, the second effect dominates and the cross section is expected to 
be larger than the conventional on-shell Born approximation. At subasymptotic energies 



the first effect is important and one overestimates the cross section. The Monte Carlo 
simulation based on the coherent branching algorithm || predicts that the b-quark lepto- 
production cross section at Hera energy is lower than the one obtained by a conventional 
one-loop calculation [^]. 

3) Coherent branching for x — > and x — > 1. 

It has been recently shown ||] that in the small x region one can resum the leading contribu- 
tions of gluon emission by a branching algorithm which has the following two characteristics. 

(i) Phase space for the branching. Destructive interference among soft gluons depletes 
the emission phase space and one finds that both for large and small x the emission takes 
place in the angular ordered region. Denoting by q ti and Zi the angle, the transverse 
momentum with respect to the incoming hadron and the exchanged energy fraction of the 
i emitted gluon we have, for small z iy 

{9 i+1 > 9i} ~ {qti+i > Ziq ti } . (4) 

For small Zi this phase space is larger than the conventional one corresponding to transverse 
momentum ordering and leading to the one-loop anomalous dimension. 

(ii) Non-Sudakov form factor. In the region x — > some of the gluons have Zi — > 0. 
The corresponding virtual corrections contain In Zj-singular contributions, which factorize 
and exponentiate to give the following non-Sudakov form factor 

\ is (zi, q ti , k ti ) == exp 

where k ti is exchanged transverse momentum resulting after the emission of gluon g^, i.e. is 
the total transverse momentum of the system formed by all partons emitted within a cone 
of aperture 9,i. This form factor becomes negligible for finite x but for small x is important 
and has to be considered together with the usual Sudakov form factor. 

The non-Sudakov form factor has the effect of screening the 1/zi singularity of the gluon 
splitting function 

A ns (zi,fti, k ti ) jzi -> Zi->0. (6) 

The /^-dependence in A ns makes the branching non local, i.e. dependent on the develop- 
ment of part of the emission process. Because of this non local /^-dependence, the new 
branching does not leads, for x — > 0, to the Altarelli-Parisi equation for the structure func- 
tion. One obtains instead the Lipatov equation for the structure function which gives the 
"Lipatov" anomalous dimension. 

Angular ordering is also the phase space constraints obtained from coherence of soft 
radiation in the large x region ||. Therefore it is possible to construct a unified coherent 
branching algorithm valid both for x — > and x — > 1, which for x — > takes into account 
the mentioned results to all-loops, and for x finite takes into account all the leading contri- 
butions and the next-to-leading corrections important in the large x region. The important 
difference between the time-like and the space-like branching is the presence in the latter 
of the non-Sudakov form factor. 



n \z t / \z t q t i 



(5) 



4) Unified equation for the structure function. 



From the coherent branching algorithm valid both for small and large x, one obtains 11 
a unified equation for the generalized structure function J-(x, Q, k t ) giving the probability 
(per unit of lnx) of finding a gluon at longitudinal momentum fraction x and transverse 
momentum k t in a hard process at the scale Q. This hard scale gives the maximum available 
angle for the branching. Taking into account only the contributions which are singular for 
x — > and 1, for the energy moment distribution one obtains 

Fn(-i**-— Q)- (7) 

z z 

In the Kernel the emission of a gluon with energy fraction z and transverse momentum 
q t = (1 — z)j z Q is factorized. In the integrand we have the distribution before the emission 
of this gluon with the total transverse momentum k t — q t and at the scale Q/z which is 
given by angular ordering. The integration over the azimuthal direction of q t is understood. 
The term ()+ is the usual 1/(1 — z) singularity for the soft gluon emission. Its regularization 
corresponds to the usual Sudakov form factor. For iV > 1 we have z = 0(1) thus both the 
rescaling of Q in F N and the non-Sudakov form factor can be neglected and the equation 
becomes the usual light cone expansion evolution equation. The regular z(l — z) — 2 term 
in the gluon splitting function is not included in this approximation. For iV — ► 1 the 
rescaling of Q and the screening of the 1/z singularity by the non-Sudakov form factor 
are important. The equation has not any more the structure of an evolution equation and 
becomes equivalent to the Lipatov equation and gives the anomalous dimension in Eq. (jH). 

5) Monte Carlo simulation. 

From the QCD coherent branching algorithm valid in all regions of x, we have constructed 
|,| a Monte Carlo simulation program which has been applied || to the study of heavy 
flavour leptoproduction at Hera and higher energies. We have compared the results of 
this program, which for small x reproduces the mentioned "all-loop" contributions, with 
the one of the conventional "one-loop" branching, which for small x reproduces only the 
first contribution of the anomalous dimension in ([D). For large x the two branchings are 
equivalent. In the following I give a short summary of the main results of Ref. |J. 

In the study of heavy flavour leptoproduction, the most important differences between 
the improved all-loop branching and the conventional one-loop branching are seen in the 
final state gluon distributions. These differences arise from the additional phase space 
available for primary gluon emission in the all-loop evolution: the region of disordered 
transverse momenta is forbidden in the one-loop evolution, whilst in the new treatment 
it is allowed, although suppressed at very small momentum fractions by the non-Sudakov 
form factor. Therefore the number of emitted gluons is enhanced, especially at small x 
and large angles, i.e. in the low-rapidity region. At present, these differences are small 
compared with uncertainties due to our lack of knowledge of the input gluon distribution. 
This underlines the importance of determining the gluon structure function experimentally 
down to the lowest possible values of x. 

One of the most important feature of the new formulation is the suppression of large 
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energy and rapidity gaps, and large neighboring pair masses, in the distributions of pri- 
mary emitted gluons. When the full colour structure of final-state branching is taken into 
account, this will have the effect of suppressing the production of high- mass colour-singlet 
combinations of partons. This in turn will produce a more local preconfinement of colour 



11| and permit a more direct connection between the perturbative parton shower and 



the observed hadron distributions. 

The inclusive distributions generated by one-loop and all-loop evolution are rather sim- 
ilar. This is mostly due to a cancellation of leading higher-order corrections in inclusive 
observables, reflected in the absence of leading singularities of the "Lipatov" anomalous 
dimension in second and third order and to the presence, for running as, of the cutoff in 
the exchanged transverse momenta. 

The effects of the leading singularities of the coefficient function, represented by the 
off-shell photon-gluon subprocess cross section, are visible in the transverse momentum 
distribution of the interacting gluon. However, the asymptotically dominant effect of the 
high-fc t tail is masked, even at very high energy, by the lower value of the off-shell cross 
section at /c^2M. Although the dynamical simplifications prevent firm predictions of the 
resulting cross sections, these results are confirmed qualitatively by recent analytical studies 
of sub-asymptotic effects Wl . 



The program || used to generate these Monte Carlo results is available from the au- 
thors. However, it should be emphasized that the program is not an event generator: it 
produces weighted events only and the weight distribution is broad. It also makes nu- 
merous kinematical and dynamical simplifications. Therefore it not suitable for detailed 
quantitative phenomenology, but is intended rather as a theoretical tool to permit a first 
look at exclusive small-x phenomena. 
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